On the determination of the last stable orbit for circular general relativistic binaries 

at the third post-Newtonian approximation 
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We discuss the analytical determination of the location of the Last Stable Orbit (LSO) in circular 
general relativistic orbits of two point masses. We use several different "resummation methods" 
(including new ones) based on the consideration of gauge-invariant functions, and compare the results 
they give at the third post-Newtonian (3PN) approximation of general relativity. Our treatment is 
based on the 3PN Hamiltonian of Jaranowski and Schafer. One of the new methods we introduce 
is based on the consideration of the (invariant) function linking the angular momentum and the 
angular frequency. We also generalize the "effective one-body" approach of Buonanno and Damour 
by introducing a non- minimal (i.e. "non-geodesic" ) effective dynamics at the 3PN level. We find that 
the location of the LSO sensitively depends on the (currently unknown) value of the dimensionless 
quantity ecstatic which parametrizes a certain regularization ambiguity of the 3PN dynamics. We 
find, however, that all the analytical methods we use numerically agree between themselves if the 
value of this parameter is astatic — —9. This suggests that the correct value of aj sta tic is near —9 (the 
precise value u-Ctatic = — 4f + fj 71 " 2 — —9.3439 . . . seems to play a special role). If this is the case, we 
. then show how to further improve the analytical determination of various LSO quantities by using 

a "Shanks" transformation to accelerate the convergence of the successive (already resummed) PN 
C^) , estimates. 
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""^J 1 I. INTRODUCTION 

cr 

\ The study of the late dynamical evolution of binaries made of compact objects (neutron stars or black holes) 
is important because such systems are the most promising candidate sources for interferometric gravitational-wave 
detectors such as LIGO and VIRGO. In particular, the global structure of the gravitational waveform emitted by such 
• rH , a binary sensitively depends on the frequency at which the system's orbital evolution changes from a gravitational- 
/\ radiation-driven inspiral phase to a plunge phase followed by coalescence • 

In the test-mass limit (// <C M) the orbital dynamics is that of a test particle (of mass fi) moving in a Schwarzschild 
background (of mass M). A very important qualitative feature of circular orbits in such a background is the existence of 
a Last Stable Orbit (LSO) located at the (area) radius i?LSO = QGM/c 2 . When considering the effect of gravitational- 
radiation reaction, one expects the test-particle motion to change abruptly near i?LSO from a slow inspiral to a fast 
plunge. By analogy, one believes that the motion of a compact binary made of comparable masses, mj and mi , will 
exhibit a similar transition from inspiral to plunge, with the location of the transition being mainly determined by 
the existence of a Last Stable Orbit in the conservative part of the two-body Hamiltonian. 

Several authors have tried to estimate the location of the LSO in comparable-mass (compact) binaries. An early 
analytical estimate was made by Clark and Eardley || using the first post-Newtonian (1PN) Hamiltonian. Some 
authors tried to use initial value formalisms to locate the LSO. However, the initial value approaches used in 

these works assume a conformally flat metric, and therefore do not correctly incorporate the well known second post- 
Newtonian (2PN) dynamics . In this work we shall take the view that a correct incorporation of all 2PN effects 
is a necessary (if maybe not sufficient) prerequisite for an accurate determination of the location of the LSO. Previous 
treatments that used the full 2PN dynamics to try to analytically determine the location of the LSO include Refs. 
p0| ^2|,p]Jlq]. The aim of the present work is to extend these PN-based analytic determinations of the LSO to the third 
post-Newtonian (3PN) level. The conservative part of the 3PN Hamiltonian for two point masses has been obtained 
in 1998 by Jaranowski and Schafer |l4| ], though with some remaining ambiguity due to the need to regularize the 
divergent integrals entailed by the use of point- like sources. The determination of the 3PN dynamics has been recently 
completed by deriving the Hamiltonian in a non-mass-centered frame, and by fixing a certain momentum-dependent 
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regularization ambiguity [ pJ5| . Recently, we have extracted from the 3PN Hamiltonian of Ref. [[Uj all its dynamical 
invariants, i.e. all the functions linking dynamical quantities which do not depend on the choice of coordinates in 
phase-space [H. In this paper, we shall use several of these 3PN invariants to determine the location of the LSO. 
Some of the methods of LSO determination that we shall use below generalize previous works PMLi, but others are 
new. 



Before embarking on the discussion of the methods we shall use here to extract the LSO from PN expansions, let 
us stress that the basic theme underlying our endeavours is the following: Our problem is to extract some semi-non- 
perturbative information from (badly convergent) perturbation expansions. We shall do that by using several types 
of "resummation methods" . The basic idea of all resummation methods is simply the following: to complete the 
information contained in the first few terms of a perturbative expansion f{z) — cq + c\ z + • • • + c n z n + 0{z n+l ) 
by injecting some non-perturbative information about the global behaviour (if possible in the complex plane) of the 
exact function f(z). The amount of global information one has about the function f(z) determines the best type 
of resummation method to use. For instance, if the only information at our disposal is that the function f(z) is 
(probably) meromorphic in the complex z-plane, then the best, all-purpose method is to use Pade approximants. If 
we knew more about the location of the singularities of f{z) in the complex plane one might contemplate to use other 
methods (e.g. change of independent variable, Borel transform, . . .). 

In this paper, we shall use two distinct classes of methods for extracting the (invariant) location of the LSO from the 
knowledge of the PN-expanded dynamics. The first class of methods (which was introduced in Ref. [Q) is discussed 
in the present Section. The second class will be discussed in the next Section. The first class of methods is based on 
the combination of three ideas: 

(i) to work only with invariant functions; 

(ii) to make a maximal use of the known, exact functional form of invariant functions in the test-mass limit 
(v = mim,2 / (mi + m^) 2 — > 0) and to assume some structural stability when the parameter v is turned on; 

(iii) to use Pade approximants to represent the invariant functions which are (because of (ii)) expected to be 
meromorphic functions of their argument. 

This method was applied in Ref. |l| to the two invariant functions which play an essential role in the gravitational- 
damping-driven inspiral of a binary system: the binding energy E(x) of a circular orbit, and the gravitational-wave 
flux F{x) emitted by a circular orbit, both being considered as functions of the (dimensionless) invariant parameter 



where M = m\ + mi denotes the total mass of the binary, and oj the orbital angular frequency along a circular orbit. 

Before introducing some variations on this method, let us motivate the interest of using the three ideas (i)-(iii) 
above. First, we recall that the PN expansions of non- invariant, i.e. gauge-dependent, functions can have (and do 
have, in some gauges) worse convergence properties than the PN expansions of invariant functions. For instance, the 
PN expansion of the gravitational- wave flux, from a test mass in circular orbit, say Ftm, considered as a function 
of the harmonic-coordinate parameter 7 = GM/(c 2 rh a rmonic), is such that the 1PN "correction" to the leading 
"quadrupole" result becomes fractionally larger than 100% (while being negative!) for a radius Tharmonic larger than 
the LSO (which is at ^lso = 5GM/c 2 in harmonic coordinates). More precisely, if one formally writes down the 
expansion of -Ftm (7), in powers of 7, near the LSO (i.e. for 7 near 1/5) one gets a series numerically proportional 
to: 1 - 1.74 (5 7 ) + 1.12 (5 7 ) 3 / 2 + 1.29 (5 7 ) 2 + • • •, i.e. a series whose first terms do not exhibit any convergence near 
the LSO. By contrast, the flux Ftm expanded in terms of the invariant parameter x, Eq. ( |2.1[ ), (xlso = 1/6 in the 
test-mass limit) has a more reasonable expansion proportional to 1 - 0.619 (6x) + 0.855 (<oxf /2 ~ 0.137 (6x) 2 + ■■■. 
Still, it is clear that one needs some resummation technique for summing a series as slowly converging as Ftm {%) . It 
was shown in great detail in Ref. JD], by making use both of the known analytical results on high-order (5.5PN) terms 
in the post-Newtonian expansion of the test-mass flux function Ftm {%) [P-7[ , and of the existence of a pole- like blow 
up of Ftm (a?) at x — a^ght ring = 1/3, that one could considerably speed up the convergence of the straightforward 
(Taylor-like) PN series, by replacing it by a suitably defined sequence of Pade approximants (see, notably, Fig. 3 in 



To further illustrate the idea (ii), and the need for an acceleration of convergence, let us recall the treatment of the 
energy function E{x) introduced in [Q. In this paper we follow p| in defining the dimensionless energy function 



II. METHODS DIRECTLY BASED ON DYNAMICAL INVARIANTS 
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where £ R is the total ( "relativistic" ) energy of the binary system (including the rest-mass contribution), and where 

», . _ m 1 m 2 _ fJ, m 1 m 2 . , 

M = mi+m 2 , \i= ■ , v= — = - ■ -3. (2.3) 

mi +m2 M (mi +7712) 

For notational simplicity we shall henceforth use units such that c = 1 . Note that the energy function used in (!]] was 
e dis = ( £ r _ m)/M = vE hcic . 

Because the argument x is, from its definition (pj|), of formal order 0(c~ 2 ), the knowledge of the dynamics at, say, 
the nth post-Newtonian (nPN) order entails the knowledge of the expansion of the ratio E(x)/x up to x n : 

E{x\v) = -iar [1 + E x (v) x + E 2 (v) x 2 + ■ ■ ■ + E n {y) x n + 0(x n+1 )] . (2.4) 

[The term —\x corresponds to the Newtonian binding energy £ NR = — i [iv 2 . Th en the term Ei(v) x, for instance, 
represents the fractional 1PN effects, etc.] The symmetric mass ratio v, Eq. ( |2.3| ), enters the expansion coefficients 
E n (y) as a parameter. At present, only E\(y), E 2 {y) and E^iy) are known (with some ambiguity for the 3PN 
coefficient E^{v)). They were written down in [jl6| and will be repeated below. On the other hand, in the test- 
mass limit v — > 0, the coefficients E n (0) are known (in principle) for any n. What is more, the exact expression of 
E[x\ v = 0) is known: 

E(x;u = 0)= (2.5) 
VI — 3x 

It is known (see, e.g., [Q) that the location of the LSO corresponds just to the minimum of the function E{x) 
((dE(x) / dx) XljSO = 0). Therefore, it would seem that the most straightforward way of locating the LSO is to consider 
the successive "Taylor approximants" of E(x), say Er n (x) (defined for each integer n as the R.H.S. of Eq. (|2.4|) , 
without the O (x n+1 ) error term), and to solve the equations dEx n (x)/dx = 0. Let us see what this gives in the 
test-mass limit where the Taylor expansion of Eq. @.5|) is known: 

1 (, 3 27 o 675 , 3969 4 \ 

E(x\ u = 0) = — x 1 - - x x 2 x 3 x 4 . (2.6) 

1 ' ' 2 V 4 8 64 128 / v ' 

The successive "Taylor" estimates of xlso> or better ^lso/^lso* = ^^LSO' are f° unc i to be (in the test-mass case): 

6x Tl =4; 6 x T2 = 1.49284; 6 x n = 1.17565 ; 6 x Tl = 1.07680 . (2.7) 

[To avoid confusion note that we use here the convention that "TV corresponds to the nPN approximation, i.e. a 
(i>/c) 2n -accurate result, while in Ref. 0] 'T n ' referred to (u/c)"-accuracy, i.e. to the -|PN approximation.] From Eq. 
(2.1) the corresponding values of the orbital frequency at the LSO, scaled to the exact value, uj = uj/u; cx£ict = (6a;) 3 / 2 , 



read 

Q Tl = 8 ; Q T2 = 1.82398 ; cD T:5 = 1.27472 ; uj Ti = 1.11739 . (2.8) 

As the value of the frequency at the LSO is the most important observable one wishes to know (for data analysis 
purposes), one should reject any method which does not have the prospect of determining it to, say, better than about 
10%. The test-mass results^ ( p.8[ ) suggest that, if the dynamics is known only up to the 3PN level, straightforward 
Taylor approximants of the energy function E{x) do not converge fast enough to determine satisfactorily the location 
of the LSO. [The 4PN level might barely suffice, but it seems anyway excluded that one will be able to analytically 
derive the 4PN dynamics. 



This preliminary discussion motivates the necessity of boosting up t he co nvergence of the series ( p.4| ) or (2.6). This 
is here that the ideas (ii) and (iii) enter. The exact test-mass result ( |2.5D suggests (under the assumption that the 
v 7^ case represents a structurally stable deformation of the v = limit) that the function E(x; v), for v ^ 0, will 
have a branch-cut singularity (at some point xq = ^+0{v)) in the complex x plane. It is, a priori, much better to work 
with functions which are meromorphic in the complex plane, because we can then make use of Pade approximants, 



1 The logic here is to use the known convergence properties of the v — limit to estimate the convergence when v 7^ 0. It is, 
indeed, unlikely that turning on v will drastically improve the convergence properties. 
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which are efficient tools for accurately representing meromorphic functions. This suggests to work with some (to be 
defined) new invar iant energy function that is known to be meromorphic in the test-mass limit. 

Looking at Eq. (2.5) one is tempted to consider the square of the function 1 + E. However, Ref. Jl| remarked that 



the usual test-mass limit definition of this function, namely (m-2 <C m\ being the test-mass orbiting mi) 

pR 

l + £!=_2_ (2.9) 

m 2 

with 

£ t R ot = m x + S* + 0(m\) =m x - + (m\) , (2.10) 

mi 

looks unnaturally asymmetric with respect to the labels 1 and 2. [Essentially, this asymmetry comes from the fact, 
hidden in the symmetric definition ( [2.2] ) , that £ 2 represents the sum of the rest-mass energy of m.2 alone and of the 
binding energy of the binary system.] Ref. jjj therefore suggested to introduce the more symmetric function 

2 2 /r-R\2 2 2 

_ s-mj- mj _ {£ n Y -m,j- m,j 
2mim 2 2miW2 

of the Mandelstam invariant s = — (Pi + P%) 2 = (£ R ) 2 ■ Indeed, in the test-mass limit 

via) = -^L ^ £ ± = l + E = ±==*== . (2.12) 

mi m 2 m 2 v 1 - 3a; 

It is then natural to define the function e(x) (or rather 1 + e(x)), by setting 

l + ei X )^(v( S ( X ))) 2 =( i£R) l- ml - m2 A . (2.13) 

\ 2 mi m.2 / 

Then, 1 + e(x), being equal to (1 — 2a;) 2 / (1 — 3a;) in the test-mass limit, i.e. (after subtraction of the trivial constant 
1) 

1 - 4a; 

e(x;u = 0) = -x — — -, (2.14) 

1 — 3a; 

one expects the function e(x\ v) to be meromorphic in x when v 7^ (with a pole located at some a: po i = \ + 0(v)). 
This finally leads to the following "P-approximant" -improved method for locating the LSO: starting from the 



Taylor approximants (T-approximants) of the original E(x) function, Eq. (2.4), compute first the corresponding 
T-approximants of the new e(x) function, say 

e(x) = -x [l + ei(i/) x + e 2 {v) x 2 + e z {v) a; 3 + C(a; 4 )] . (2.15) 

Then construct a sequence of Pades of the Taylor-expanded e(x): 

epjx) = P e k [T n [e(x)]] , (2.16) 

where k + £ = n. Here k and £ are the degrees of the polynomials Nk(x) and Di(x) entering the Pade Pi{x) = 
Nk(x)/ Di(x). It is known that, generically, the Pade improvements are best when one is near the "diagonal", i.e. 
when \k — £\ is "small" compared to k and £. When dealing with a function f(x) that is expected to have a pole at 
some xo 7^ 0, one imposes the constraint £ > 0. At 1PN order this uniquely fixes the values of k and £, namely k = 
and I = 1. At 2PN order the Pade closest to the diagonal is that with k = 1 and 1 = 1. At 3PN order there are two 
possible Pades near the diagonal, namely k — 1 and I = 2, or k = 2 and £ = 1. In this work we shall use the latter 
one (Pi) because we found it to be more robust under variations of the coefficients of the Taylor expansion of the 
Padeed function. (One aspect of this robustness is that the existence of a real pole in this Pade is always ensured, 
while this is not the case for the other 3PN possibility: Pj.) 

Note also that Pades are originally defined only for series of the regular type a(x) = Co + ci x + ■ ■ ■ with Co 7^ 0. 
When dealing with a function of the type f p (x) = x p o~(x) with some relative integer p, we shall, by convention, define 
any Pade of f p (x) as being (k + I = n) 
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P?[T n \f p (x)]] ee a*P?[T n [x-* f p (x)}] . (2.17) 

The "ep n -estimate" of the location of the LSO is then defined as the value xp n at which ep n (x) reaches a minimum. 
[It is easily seen that e(x) follows the variations of E(x), and in particular that it reaches a minimum at the same 
place as E(x).} 

We spent some time explaining in detail on one example our general Pade-improved-test-mass-limit-motivated 
approach because we are going to extend it to several other invariant functions. In this (and the next) section, we 
present our various methodologies. The 3PN results obtained by them will be presented in a later section. 

The introduction of the function e(x) has two defects. First, it is not unique because we do not know for sure 



whether the function (<p(s)) 2 , Eq. (2.13|), is a "better" invariant than (1 + E) 2 = (1 + (y/s — M)/fi) 2 . Second, the 



Padeing of e(x) starts giving meaningful results only at the 2PN level. Indeed the 1PN expansion of e(x), in the 
test-mass limit, e(x; v = 0) = — x (1 — x + 0(x 2 )), yields a Pade ep 1 (x; v = 0) = — x/ (1 + x) which contains no pole 
on the positive real axis, and which formally predicts an LSO (minimum of ep 1 (x)) located at x = +oo. 

In this paper, we propose to consider another invariant function, which is more uniquely defined, and which gives 
sensible results already at the 1PN level. Let us consider the reduced angular momentum 

j = ^ J_ ^ (2.18) 

/iGM Gm\m2 ' 

where J denotes the total angular momentum of the syste m. I n the test-mass limit the invariant function giving the 



(dimensionless) quantity j in terms of the quantity x, Eq. (2.1), reads 



j(x;v = 0)= 1 (2.19) 
V x(l — dx) 

This motivates the consideration of the squared (reduced) angular momentum j 2 (x) which is expected, when v ^ 0, 
to be a meromorphic function of x, with a pole at the "light ring" x po io = \ + 0(y). Therefore we propose to work 
with the Padeed form of j 2 : 

j Pn (x;v) = P£[T n \j 2 (x;v)]}, (2.20) 

with k + £ — n, and the choice of £ > discussed above. (We use in Eq. ( [2.20| ) the convention ( [2.171 ), i- e - the factor x _1 
in j 2 (x) is factored before taking a Pade.) Note that, in the test-mass limit, if we knew only the 1PN approximation 
to the function j 2 (x), i.e. j 2 (x; v = 0) = x _1 (l + 3x + C(x 2 )), the procedure ( [2.20| ) would reconstruct the exact result: 
j 2 i (x;^ = 0) = P 1 °[T 1 [i 2 (x^ = 0)]] = [x(l-3x)]- 1 . 

It is important to note that the perturbative information contained in the PN expansion of the function j(x) (or 
equivalently j 2 (x)) is totally equivalent (at any PN accuracy) to the information contained in, cither the original 



energy function E(x), or the new one e(x), Eq. (2.13). Indeed, the generic Hamiltonian equation 9i = = dH/dli 
in action-angle variables (Ii,9i) yields 



^circular 



d Il-^ d _E (221) 

dj " GM dj ' 1 ' 



which implies the identity 



dE(x) _ 3/2 dj(x) 



x 6lz ^- L . (2.22) 



dx dx 



The identity ( 2.22 ) proves the assertion just made about the identical information content in E(x) and j(x). It 
also proves several interesting facts. First, the location of the minimum of (the exact) j(x) coincides with that of the 
minimum of (the exact) E(x) (both of them equivalently defining the LSO). Second, the existence of a branch cut 
singularity oc (xo — x) -1 / 2 in either j(x) or E(x) necessarily implies the presence of a similar singularity oc (xo — x) -1 ' 2 
(at the same location xo) in the other function (E(x) or j(x), respectively). This can be viewed as a confirmation of 
our generic assumption of structural stability. However, this argument also shows the ambiguities present when trying 
to work with the energy function. Indeed, if we assume that, near xq, j(x) can be expanded as ip(x)(xo — x) -1 / 2 , 



where <p(x) is a smooth function, one finds from Eq. ( [2.22 ) that E(x) = i(}(x)(xo — x) x l 2 + c, with some unknown 



constant c. The lack of knowledge of the constant c (which a priori depends on the parameter v) implies that we do 
not know which function [E — c(v)) we would square to get a meromorphic function with a simple pole. The same 
reasoning shows another defect of the proposal to consider the (new) function e(x). Indeed, when E(x) — > oo the 
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leading term in e(cc), Eq. ( pl3| ), will be e{x) ~ jV 2 E 4 (x), which will have a double pole oc v 2 (x§ — x) 2 , ii j 2 {x) has 



a simple pole oc [xq — x)~ x . For all these reasons, we consider that the "j-method", Eq. ( 2.20| ), appears as the best 
way of locating the LSO, within the class of methods dealt with in this section. 

To conclude this section, let us, however, mention another invariant function one might wish to consider. This 
function is the fourth power of the dimensionless periastron parameter 1 + k = $/(2tt) considered as a function of the 
reduced angular momentum. Indeed, in the test-mass limit, and for circular orbits, one knows that |||l8],[l6| 

(l + ^) 4 =(l-^) • (2-23) 

We recall that j 2 = 12 is the location of the LSO. Therefore, if we define K = (1 + fc) 4 = ($/(27r)) 4 and y = 1/j 2 , we 
might consider 

Kp n {y;v)=P^[T n [K{y-v)}], (2.24) 

with k + I = n, and our canonical choice of t > 0. Then, we can take the pole of Kp n (y; v) as estimate of the value 
of 1/j 2 at the LSO when v =^ 0. We consider, however, that the j-method (or any energy method for that matter) 
has a better chance of accurately locating the LSO because it incorporates not only the information that something 
special (a minimum) takes place at the LSO, but also the information that the location of this minimum is a corollary 
of the presence of a blow up (j 2 — > oo, E — > oo) at a point, x po ic, further away on the real axis. Therefore, even if the 
location of x po \ c is not known very accurately, one can hope that cclso will be more robustly determined. 



III. EFFECTIVE ONE-BODY METHOD 



In this section, we shall turn to a rather different method, though one which also makes use of the three ideas 
(i)-(iii) listed at the beginning of the previous section. This new method incorporates a fourth idea, which has been 
recently put forward by Buonanno and Damour fl3j | . This fourth idea consists in mapping (through the use of invariant 
functions) the real two-body problem we are interested in (two masses mi, mi orbiting around each other) onto an 
"effective one-body problem" (one mass mo moving in some background metric, ff^ cctlvc (^ 7 ))- At the 2PN level, 
Ref. HU has shown the possibility of mapping the real two-body problem onto geodesic motion in some spherically 
symmetric metric g^ ctlvc (x 1 ; v). It was found that, when v ^ 0, ff^ ectIVC (a; 7 ; v) is a smooth deformation of the 
Schwarzschild metric. The "mapping rules" between the two problems were motivated by quantum considerations: 

(i) the adiabatic invariants Ii — § Pi dqi (which are quantized in units of fi) were identified in the two problems, and 

(ii) the energies are mapped through a function ^effective = / [frcai] which is, a priori, arbitrary, and which is 
determined in the process of matching the two problems. 

In other words, the idea is to determine a metric g^a ctlvc such that the "energy levels" £ effective [^i] (i.e. the Hamil- 
tonian expressed in action variables, or "Delaunay Hamiltonian" ) of the bound states of a particle moving in g^ ectlve , 
are in one-to-one correspondence with the two-body bound states: 

^effective [Ii] = f [£rcal [Ii]] ■ (3.1) 

The unknowns of the problem are the numerical coefficients entering a spherically symmetric metric 

ds 2 cS = -A(R eS ) dt 2 cff + dR l s + R 2 S {d9 2 cB + S in 2 9 cB d^ ff ) , (3.2) 

A.{HeS ) 



namely, 



GM [GM Q \ 2 {GM \ 3 /GAV 4 
— = h a? — — + as — =r— + cu — — 



A (j R) = l + ai _^ + a2 ^_^j +a3 ^_^j +a4 ^_^j +..., (3 .3a) 

^) = i+^V + HVJ +d3 {~ir) + "' ! (3 - 3b) 



and the coefficients entering the energy-map / (here written for the "non-relativistic" energies £^ = £ G ff — ttiq, 

NR. 
eal 



£™i = freal " M): 



G 



cNR cNR 

c cS _ c rcal 



TO 



1 + ai 



cNR 
c rcal 



a 2 



cNR. \ 2 / cNR 

c rcal \ i / c rcal 



a 3 



(3.4) 



As discussed in Ref. fll3|j it is natural to require that the effective mass too be exactly equal to the usual non- 
relativistic effective mass /i = m\m%/{m\ + to 2 ). One can also (by convention) choose the mass Mo entering the 
effective metric to be Mq = M = toi + to 2 . With these choices the Newtonian limit tells us that the first coefficient 
in A(R) = —goo is simply a\ = —2. The 1PN level then involves the coefficients a 2 , d\, and a±, while the 2PN 
and 3PN levels involve (a.3, <i 2 , a 2 ) and (04, d^, a^), respectively. In other words, at each PN level, we only have three 
arbitrary coefficients to play with. This seems to be quite a small number of degrees of freedom, compared to the many 
possible coefficients which can enter the PN-cxpansion of the Delaunay Hamiltonian. [The Delaunay Hamiltonian 
was determined at the 2PN level in Ref. H, and at the 3PN level in Ref. pq|.] In order to clarify the number of 
independent equations to be satisfied when mapping the real problem onto the effective one, let us consider a generic^] 
PN-expanded Hamiltonian, with the symbolic structure 



-p 



2n-2 



(npf 



1 



M 2 1 + Z2 



2(n-l) 



(np) 



2(n-l) 



1 
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n+l 



(3.5) 



Here, we consider the reduced Hamiltonian 7J NR = iJ NR //i, in the center-of-mass frame, as a function of the reduced 
canonical variables p = Pi/V = — P 2 /V, q = (xi — x 2 )/(GM) ((np) denotes n • p with n = q/g). See Ref. [[l6| for 
details. Note that we follow here [13| in denoting by q, p the original (ADM-like) coordinates. We have suppressed 



all coefficients in Eq. (^5) to display the st ruc ture. What is important for our present purpose is the total number of 
coefficients in the nPN- level Hamiltonian (pTa). This is easily checked to be: 



C H (n) 



(n + l)(n + 2) 



+ 1 



(3.6) 



As explained in Ref. one way (and the only explicit one) to map the real Hamiltonian ( |3.5[ ) onto an effective 
Hamiltonian, while keeping the action variables invariant, is to apply a canonical transformation, with generating 
function G(q,p') — q l p[ + G(q,p'). The most generic generating function that we need to consider has the symbolic 
structure (at the nPN level) 



GnPN(q,p) = (q p) \p 



2(n-l) + . . . + ( np )2(n-l) 



1 



(3.7) 



Correspondingly to the pure p-dependence of the leading kinetic term in fl3.q ) we have written here the leading 
term in G„pn as cx (qp)p 2n - [It is easily shown that any term of the form (qp) p 2 ("~ fc ) (np ) 2k must have a vanishing 
coefficient.] The number of arbitrary coefficients in the nPN-level generating function (|3.7|) is easily seen to be 



C G (n) =C H (n-l) 



n(n + 1) 



1 . 



(3.8) 



Finally, the difference A(n) between the number of equations to satisfy, and the number of unknowns (including the 
3 basic parameters (a„ +1 , d n , a n ) appearing in the effective metric and the energy- map) reads, at the nPN level 



A(n) = C H (n) - C G (n) 



(3.9) 



In particular: A(l) = —1, which means that requiring a 1PN matching leaves one degree of freedom unrestricted. 
[This freedom was used in [fl3f to impose the natural condition d\ = 0, i.e. that the linearized effective metric coincides 
with Schwarzschild.] Then A(2) = 0, which means that there will (barring any degeneracy) be a unique solution at 
the 2PN level. [This was indeed the result of [l3|.] But A(3) = +1, which means that, at the 3PN level, there is one 
more equation to satisfy than the number of free parameters. Then the situation would get worse and worse at higher 
PN levels. 

By exp licitly performing the matching between the canonically-transformed Hamiltonian and (modulo the energy 



map (3.4)) the effective Hamiltonian of a point particle moving in some gf^,, we have established (details will be 



We use the information that the leading kinetic terms in a PN Hamiltonian are given by the expansion of the free Hamiltonian 
\/pi + m\ + \/p2 + m%, so that, at the nPN level, they are cx p 2 ( 1+n ' without dependence on n = q/q. 



7 



given below) that, if we follow Ref. [T^] in imposing the natural condition d% = at the 1PN level, there are, indeed, 
Ch(3) = 11 linearly independent equations, for Cg(3) + 3 = 10 unknowns, to be satisfied at the 3PN level. [No miracle 
occurred!] 

At face value, this looks like bad news for the idea of the "effective one-body approach" . However, we think that 
there are acceptable ways to rescue this approach. A first cure would be to take advantage of the fact that the total 
number of equations at the three first PN levels (1, 2, and 3) is exactly equal to the number of unknowns (in other 
words A(l) + A(2) + A(3) = —1 + 0+1 = 0). Therefore if we relax the (not really nece ssary) constraint d\ = 0, 
there will be a unique 3PN-accurate effective metric g^ v (and a unique energy mapping (3.4)) satisfying the necessary 
constraints. We have verified that this is true and, for completeness, we give this unique solution in Appendix A. 
But, we do not want to take this solution too seriously for the following reasons: (i) it does not look natural to have 
to wait to know the 3PN Hamiltonian to determine the 1PN and 2PN effective metrics; (ii) this solution looks more 
complicated than the 3PN Hamiltonian itself; and (iii) this trick is not expected to be sufficient to ensure the existence 
of solutions at higher PN levels. (Indeed, A(n) = n — 2 continues to increase.) 

We propose therefore to consider a second (more radical, and simpler) way to cure the problem. Indeed, we have to 
face the fact that there is (probably) nothing deep in the effective-one-body approach. After all, it is just a somewhat 
artificial way of mapping the complicated two-body dynamics on a simpler one-body dynamics. There is no reason to 
assume that the one-body dynamics can, to all orders, be considered as equivalent to a simple geodesic motion. Let 
us recall that, in quantum mechanics^], geodesic motion means a simple Klein-Gordon Lagrangian 



fg^i ((V^) 2 + ml(p 2 ) with (Vip) 2 = g^ d a (pdp ip . 



(3.10) 



It is well-known that effective actions generally include, at higher orders, some higher-derivative terms: for instance 
of the type (□<, (p) 2 , (V</?) 4 , etc. Coming back to the classical limit, i.e. to the Hamilton- Jacobi equation (obtained 
by considering that <p(x) = exp (i S(x)/h) with h — » 0), we should correspondingly expect that, at higher orders, the 
effective one-body "Hamilton- Jacobi" equation be of the generalized form 







9cS ( X ) Va PP + A af}lS (x) p a pp p 1 ps 



(3.11) 



with p a = d S{x) / d x a . If we were to use a Lagrangian formulation, the general form (3.11) would correspond to an 



[l + A a ^s (x) u a u@ u 1 u s + • • •] with u a = dx a /ds c s, i.e. to a general (perturbative) Finsler 



action S = —m Q J ds c g 
structure. 

If we use perturbatively the lowest-order "on shell" condition (i.e. m 2 + gZ& PaV d — 0) we can (for instance) 
eliminate the presence of the energy £ c g = —po in the quartic (and higher) terms in ( |3.11 ). In other words, we can 
restrict ourselves to considering purely spatial high er-order tensors A af3 " )S (x) = A l; > M (x), etc. Dimensional analysis 
shows that the quartic terms (C(p 4 )) in Eq. (3.11) which can enter at the 3PN level must have a q-d ependence of 
the type A yfc6 (x) ~ q~ 2 . Finally, we are led to consider, at the 3PN accuracy, (after solving Eq. ( 3.11 ) with respect 
to the effective energy £ c q = —po) a generalized effective Hamiltonian of the form 



i? c R ff (q',p') = \A(q') 



( AW) 
\D(q>) 



1 ) (n' ■ p') 2 + — (zi(p' 2 ) 2 + z 2 p'2(n' • p') 2 + z 3 (n' • p') 4 ) 



(3.12) 



where the quartic (in p') terms come from the A a @ lS coupling and modify the normal "geodesic" Hamiltonian 
\J ~ -9oo ( 1 + ff off Pi Pj ) ■ I n anticipation of the need to transform (via a canonical transformation) the original coordi- 
nates (q, p) of the real (reduced) Hamiltonian into the coordinates of the effective dynamics, we have denoted the 
latter by (q', p'). 

This procedure introduces three new arbitrary degrees of freedom at the 3PN level: z\ , z 2l and 2:3. It is clear that it 



now becomes possible to map the real dynamics on the generalized effective dynamics (3.12). This becomes, in fact, 
possible in many ways. As we are primarily interested in (quasi-)circular orbits we shall find convenient to consider 
only the simple case where z\ = z 2 = 0, i.e. to use only Z3 ^ as new degree of freedom. [This degree of freedom then 
disappears in the discussion of circular orbits, which can then be considered as following essentially from a "geodesic" 
dynamics.] 

An important feature of our proposal ( [3.11 ) is that it is clearly general enough to allow for the existence of solutions 
at arbitrary PN orders. For instance, at 4PN we would have the freedom to introduce either arbitrary sextic terms, 



We recall that the origin of the effective-one-body approach lies in the quantum electrodynamics of two-charge systems, see 

isl. 
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q' 2 \p m + p' 4 (n' p') 2 +•••], or a modification of th e quartic terms: q' 3 [p' 4 + ---]. This is more freedom than is needed 
to compensate for A(4) = 4 — 2 = +2, Eq. (|3.9| ). We can also clearly arrange things so tha t circ ular orbits always 



follow from a "geodesic" dynamics. The only somewhat unsatisfactory feature of the proposal (3.11) is that we cannot 
offer any principle for determining a priori the structure of the "post geodesic" terms. In particular, it would have 
been nice to say that (as is often the case in effective actions) the additional terms are somehow generated by the 
leading term, i.e. by the effective metric. We have in mind here relations of the type: 

AafiyS = Ai R a p RfS + A2 V Q /3 R^s + A3 V Q/ 3 7< 5 R H , 

where R a p is the Ricci tensor of 3°^. However, we have checked that such "geometrical" tensors do not yield possible 
3PN corrections, but only much smaller contributions (starting at 5PN). 

We shall give in the next section the details of the computation of the coefficients entering A(q'), D(q'), the 
generating function G3PN, the energy-mapping function and Z3. For the time being, let us only stress some conceptual 
points. First, it is remarkable that the unique solution of the 3PN real <^ effective matching problem leads to the 
simple value 

a 3 = 0, (3.13) 



for the 3PN parameter entering the energy-mapping (3.4). Ref. [|13|] had already found that a-i = at 2PN, and 
ai = i//2 at 1PN. These values correspond exactly to 

= (£> (Sreal) = „ ■ (3.14) 

mo lm\mi 

We find remarkable that the simplest, symmetric function of the Mandelstam invariant s rea i = (£^ a i) 2 which reduces 
to £o/mo in the test-mass limit turns out to define the unique energy-map needed to l ink t he real 3PN dynamics to 



the effective dynamics. We interpret this as a good sign for our generalized dynamics ( 3.11 ). [By contrast, the other 
proposal of relaxing the natural constraint d\ = leads to a very complicated energy-map with ot\ ^ vjl, 012 7^ 0, 
and 0:3 7^ 0, see Appendix A.] 

Let us now consider the consequences of the effective-one-body approach for the determination of the LSO. For 
the case of circular orbits the effective-one-body approach boils down to saying that the real energy is the following 
function of the effective-one, 



£& = Wl + 2ii*LJ£ (3.15) 



(as obtained by inverting Eq. (3.14)) and that the effective energy along circular orbits is the square- root of the 



minimum value of a certain "effective radial potential" : 

y = /iRWiU, (3-16) 
where Wj(q') is obtained from (the square of) Eq. ( 3.12| ) by setting n'-p' = (and p' 2 = (n' x p') 2 + (n'-p') 2 = j 2 /q' 2 ) 

W j {q l )=A{q')(l+^+z x ^) . (3.17) 

As said above, we shall assume (as is always possible) that z\ = 0, so that the effective potential has the usual 
"Schwarzschild-like" value — <?oo((?') (1 + j 2 /q' 2 )- The value of the metric coefficient A(q') = —goo(q') is obtained, at 
the 3PN level, as a perturbative expansion in l/q' = GM/R c s: 



q <7' J cp \q' 

Note that, finally, in this approach the entire effect of the 3PN dynamics (for circular orbits) is contained in the sole 
coefficient oAv) (whose value will be discussed in the next section). 

In Ref. |13| the metric coefficient A(q') was used (at the 2PN level) as a simple truncated Taylor expansion: 
^2Pn(?') = 1 — 2/q' + 2v/q' 3 . This simple-minded approach is not adequate for dealing with the 3PN level. Indeed, 
we shall see in next section that the coefficient 04 (u) is positive and can be relatively large. In keeping with the spirit 







of the present work where we systematically try to use adequate resummation methods to improve the convergence 
of PN expansions, we shall define our effective-one-body radial potential, at the nPN accuracy (expressed in terms of 
the convenient variable u = 1/q') 



W J p »(u)=A Pn (u)(l+j 2 u 2 ) 



by using a suitable Pade approximant of Eq. ( |3 . 18 ) : 

A Pn (u) ee P k g [T n+1 [A(x 



(3.19) 



(3.20) 



with k + £ = n + 1 (because it is which corresponds to the nPN level) and, now the constraint k > (instead 

of £ > as above), because we want to factor a zero of A{u) (and no longer a pole). The Pade improvement of A(u) 
is really needed (and makes a difference) only at the 3PN level. We have found that the most robust Pades (under 
variation of the Taylor coefficients) are defined by taking k = 1 and £ = n. 

Summarizing the present method: The effective radius q' of circular orbits is obt ained as a function of the reduced 
angular momentum j by looking for the minimum of the radial potential ( 3. 19] ), where u = 1/q' and where the 
Pade-improved function A is given by Eq. (3.20), with k = 1 and £ = n. For each value of j above som e threshold 
Jmin, Wj(u) admits a (unique) minimum u*(j). From this one then determines the effective-energy ( |3.16|) , and then 
the real one ( |3. 15 ), namely 




The real circular orbital frequency corresponding to j is then given by using the identity (2.21). This yields 



GM WrealC?) 



^/TTf^U) \/i + 21/ v/w~Mj))-i 



(3.21) 



(3.22) 



Note that in this approach j = j rea i = ^effective- Finally, the LSO is invariantly defined as the minimum value of j, 
Jlso = imin, for which Wj(u) admits a local minimum. When j < j'lsOi Wj(u) has no local minimum and there 
are no (stable or unstable) circular orbits (see, e.g., Fig. 1 of Ref. [^3|). If one is only interested in locating the LSO 
(as a function of v) it suffices to look for the existence of an inflection point of Wj(u), i.e. to solve the simultaneous 
equations dWj{u)/du — and d 2 Wj(u)/du 2 = 0. 



IV. RESULTS 



Let us now give the details of the application, at the 3PN level, of the methods explained above. We follow the order 



of presentation given in the previous two sections. In what follows, we use the results of Ref. |16 for the dynamical 
invariants of the 3PN two-body Hamiltonian. We recall that the 3PN Hamiltonian derived in Ref. jlij contained two 
ambiguous parameters w s t a tic and akinetic (these ambiguities arise because of the need to regularize badly divergent 
integrals |IJ,^,^6|), and that all the dynamical invariants involve only the combination a (y) = w s t a tic v + ^kinetic v 2 

E- .... . « ^ 

Recently the 'kinetic' ambiguity parameter Wkinetic was uniquely determined |15j (see also [|21|) to be Wkinetic = 41 /24, 
so that a (y) = ^static v + % v 2 and the remaining 3PN ambiguity is embodied in the product w s tatic v. As discussed 
in the Appendix A of Jl^] one expects (both by judging from the other coefficients, and by looking at some of the 
sources of ambiguity) that w sta ti c varies in the range 

- 10 < Static < +10. (4.1) 



A. e-method 



For self-containedness let us quote the results obtained in our pre vious paper |lq| for the link b etw een the energy and 



the x-variable, Eq. (2.1). The original energy function E, Eq. (2.2), admits the PN expansion (2.4) with coefficients 



10 



where 



Ei(y) = - — (9 + v), 



E a (v) = - — (81-57v + J), 

E%{v) = - y( W i( U ) _ "static v), 



. . 405 1 / 2 6889\ 31 2 7 , 

wAv) = 1 41-7T \v -\ v 1 H iA 

w 128 64 V 6 / 64 3456 



(4.2a) 
(4.2b) 
(4.2c) 

(4.3) 



Correspondingly to this expansion, the new energy function e, Eq. (2.13), admits the expansion (2.15) with coefficients 



e\(v) 
e 2 {v) 
e 3 {v) 



1+ r 



„ 35 

3 v , 

12 ' 



10 



(w 2 (u) ~ ^static f), 



where 



1 



. . 27 1 /41 , 4309\ 103 a 

U> 2 (f) = 1 7T 1/ H V V 

K ' 10 16 V 4 15 J 120 270 



The 2PN and 3PN Pades of e(x) are given by (see Ref. @ for the 2PN case) 

l+|f- (4- + ^ 2 )x 



e P2 (x) = Pl[T 2 [e(x) 
e P3 (x) = P? [T 3 [e(x) 



35 
12 



1 - (l + \v + w 3 {y))x - (3 - f|i/ - (1 + |tj ^(4 
1 - w 3 (^)x 



where 



w 3 (j/) 



40 



" (w 2 {v) - ^static y) ■ 



(4.4a) 
(4.4b) 
(4.4c) 

(4.5) 

(4.6a) 
(4.6b) 

(4.7) 



36 - 35^ 

The corresponding x-location of the e-LSO (minimum of e{x)) can be written down analytically only at the 2PN level 



6 ^lso ( 



M 



2 - 



1 - ±L> 2 
1 16 ' 36 



(4.8) 



Note that 6 x L g 2 (i) = 1-1916, which means that the ep 2 -predicted value of the orbital frequency at the LSO differs 



from the "Schwarzschild value", GMuj[f^ = (xgg^) 3 / 2 wi t h ^lso = V6, h Y the factor 



^LSO 



^LSO 



UJ 



Schw 
LSO 



= (6x LSO ) 3 / 2 , 



(4.9) 



which is about 1.3007 in the present case. We will use Slso a s our main tracer of the "observable" location of 
the LSO. It is important to note from the start that (as emphasized in the e-method predicts (at 2PN) that 
the orbital frequency at the LSO is larger than the "Schwarzschild value". The corresponding results, at 3PN, are 
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TABLE I. Equal-mass 
binding energy Slso, Eq 
values' 



(y = 1/4) binary-system LSO parameters obtained by mea ns of different methods. 
(2.2), and the reduced angular momentum jhso, Eq. (2.18), are divided by their 1 



^lscTI — 1 ~~ 3V0 and j^fo = the dimensionless orbital frequency £Dlso is defined in Eq. (4.9) 



The reduced 
Schwarzschild 
The line 



denoted by 'BD' reports the 2PN results obtained in Ref. [13]. 







-Elso 


/l-^LSO 1 






/ -Schw 

Jlso/Jlso 








u)lso 














3PN 










3PN 










3PN 




method 


1PN 


2PN 




Wstatic 




1PN 


2PN 




^static 




1PN 


2PN 




^static 










-10 





10 






-10 





10 






-10 





10 


e-method 




-1.142 - 


1.015 


-1.369 


-1.636 




0.956 


0.994 


0.904 


0.857 




1.301 


0.998 


1.969 


2.870 


/-method 


-0.973 


-1.091 - 


1.039 


-1.322 


-1.842 


1.014 


0.970 


0.986 


0.913 


0.833 


0.960 


1.173 


1.060 


1.774 


3.915 


eff. method 


-1.007 


-1.048 - 


1.042 


-1.168 


-1.212 


1.000 


0.983 


0.985 


0.946 


0.934 


1.015 


1.075 


1.064 


1.297 


1.383 


BD, Ref. 


-1.007 


-1.050 








1.000 


0.983 








1.015 


1.079 








fc-method 












1.000 


0.980 


0.974 


0.955 


0.936 













exhibited in Table || Let us only note here that the tendency to get ljlso > 1 seems confirmed, at the 3PN level, 
rather independently of the value of the ambiguity parameter astatic- 

Once the value of £lso ( v ) is determine d (an alytically or numerically) one c an co mpute the corresponding value of 
the (real) reduced binding energy E, Eq. (§J). It is obtained by solving Eq. ( |2.13| ) in terms of £ R = M + fj,E. The 
solution is explicitly given by 



E(x) = 



\Jl + 2v (v/T+ e(x) - l) - 1 



(4.10) 



Then , knowing E(x) we can also compute the value of the reduced angular momentum j by integrating the identity 
Q2.22D, Integrating Eq. fl2.22|) by parts yields 



j(x) = -2x~ 



1/2 dE(x) 



dx 



+ 2 / dxx~ 1/2 



d 2 E{x) 
dx 2 



(4.11) 



where we have incorporated the information that j(x) ~ x^ 1 ^ 2 when x — > (i.e. in the limit of very wide circular 
orbits, described by Newtonian dynamics). By applying this result to x — xlso(^), one finally gets the value of 
Jlso iy) ■ The results so obtained by the e-method are shown in Table |[ 



B. j-method 



In this approach the basic PN expansion we consider is that of 1/ j (x). It reads (cf. Eq. (5.11) in Ref. 16|) 



1 



j 2 {x) 



1 — ^(9 + V)X + -^VX 2 - ^(wa(v) - Static v)x 3 



(4.12) 



where 



M") = ^1 4l7r 



5269 



\ 61 2 1 

v H v 2 . 

6/64 432 



(4.13) 



From Eq. (4.12) one gets 



1 + r, (9 + v) x + — (36 - v) (9 - 4i/) x 2 + ^ (w 5 [y) - Static v) x 3 
■ 36 ! 



16. 

y 



(4.14) 



where 



81 , 1/ 2 7321\ , 23 2 1 



wr(v) = 1 41 7T v H v 2 H v 

5V ' 16 64 V 6 / 64 216 



(4.15) 
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-l 



Ustatic =-10 



-1.3 
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/ 
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FIG. 1. Binary-system LSO parameter s ob tained by means of the j-method as functions of the symmetric mass-ratio 
v. The reduced binding energy £?lso, Eq. (2.2), and the reduced angular momentum j'lso, Eq. (2.18), are divided by thei r 
"Schwarzschild values": I^lscTI — 1 — \ an d jEscT = \/l2j the dimensionless orbital frequency £)lso is defined in Eq. (4.9). 
The lines shown in the plots correspond to different PN orders of approximation: 1PN (solid), 2PN (dotted), and 3PN (dashed). 



As explained above we construct the following sequence of near-diagonal Pades of j 2 (x): 

f P2 (x) = Pl[T 2 \j 2 (x)}] 
j 2 3 (x)=P 2 [T a \j 2 (x)}] 



where 



x(l - (3 + \v) x) ' 
1 



X 



(!+!„_ (3-^ +^3-)' 

1 + (3 + \v - w 6 {v))x + (9 - %v + \v 2 ~ (3 + \v) w e {v))x 2 



w 6 (u) 



192 



(36-i/)(9-4i/) 



ir(l — w%{y)x) 
(w 5 (v) - Static v). 



(4.16a) 
(4.16b) 
(4.16c) 

(4.17) 



The corresponding Xlso (now defined as the location of the minimum of j(x), or, equivalently, j 2 {x)) can be written 
down analytically at 1PN and 2PN 



1 



(4.18a) 



6 ^ ft0 - 8 ( 9 + ") 



2(9 + 



v/(36-t/)(9-4i/) 



- 1 



(4.18b) 



Note that while §x^ (y) is very slightly smaller than 1, 6x L g 2 (z^) is (like for the ep 2 estimate) larger than 1. In 

particular, 6i[g 2 (¥) — 1.1121, and the corresponding dimensionless frequency is w L g 2 (i) = 1.1728. This tendency 
to get "larger than Schwarzschild" frequency at the LSO is confirmed by the (numerical) 3PN results which are 
exhibited in Table | and Fig. |. 

Within the present j-method, once the value of xlso (y) is det ermin ed one can compute the corresponding value of 
the (real) reduced binding energy E by integrating the identity (2.22). Indeed, one can write 
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E{x) = f dxx 3 / 2 ^, (4.19) 
J dx 

where one has incorporated the boundary condition that E(x) — > when x — > 0. The results so obtained are shown 
in Table | and Fig. 0. 

C. fc-method 

For completeness , let us (though it is not among our preferred methods) mention some results obtained by using 
the fc-method, Eq. fl2.24p . The PN expansion of the function K(y), where K = (1 + fc) 4 and y = 1/j 2 , reads (cf. Eq. 
(5.27) in Ref. @) 

K(y) = 1 + 12y + 24(6 - u)y 2 + 24(w 7 (v) - ^ static u)y\ (4.20) 

where 

/41 , 128\ 1 9 , 
w 7 {y) = 72+ {^-n 2 -— y+^ 2 - (4.21) 

As explained above, we construct the following sequence of near-diagonal Pades of K(y): 

K Pl (y) ee P° [Tx[K(y)]] = (4.22a) 

Kp 2 {y) ee Pi [T 2 [K(y)}} = , (4.22b) 

1 - 2{Q-v)y 

9r r 1 + ( 12 -w s (v))y+ 12(2(6 - v) - w 8 (v))y 2 

K P M ee P 2 [T 3 [K(y)]} = * p '- BUHL, 4.22c 

1 - w$(v)y 

where 

Wa(v) EE g _ ~ ( W 7( V ) - ^static^)- (4.23) 

Then we take the poles of Kp n as estimates of the value of 1/j 2 at the LSO. The results for equal-mass binaries 
(y = 1/4) are given in Table I. 

D. Effective one-body method 

Let us explain in detail how we implemented the effective one-body method. Two methods of implementation were 
presented in Ref. jL3|. We could have used the first one by starting from the 3PN Delaunay Hamiltonian given in 
Ref. |l6| ]. However, we found finally as convenient (given the existence of good algebraic manipulators) to use the 
second method, which has the advantage of being more informative. This second method consists of writing explicitly 
the equations that have to be satisfied by the looked for generating function G(q,p') and solving them. Indeed, we 
look for a canonical transformation between the original (quasi- ADM) coordinates (q,p) of the real problem (i.e. the 
phase-space coordinates in which was obtained the orde r-reduced Hamiltonian H(q,p) in flql ), and the "effective" 



coordinates (q',p') (i.e. the coordinates used in Eq. (3.12)). The effect of the generating function G(q,p') reads 



, _ ,, , dG(q,p') , , dG(q,p') 



q H = q l + Pi=Pi+ ;T ' ( 424 ) 

ap\ aq l 

Note that, as is well known, the canonical transformation is defined only in implicit form: q' and p being given as 
functions of q and p'. But there is, in fact, no need to solve for, e.g., (q',p') as functions of (q,p). As the basic idea 
is anyway to identify the numerical value of, say, H c e(q' ,p) with the numerical value of some (to be determined) 
function f(H- [ca \(q,p)), we can do this iden tific ation by expressing both sides in terms of any set of common variables, 
say q and p' . Finally we write (using Eq. ( [j.4| )) 
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H* s ( q >( q ,p'),p')} = \l + H™( q ,p( q ,p')) 



nrV 

real J 



(4.25) 



in which the L.H.S. is given by the R.H.S. of Eq. Q3.12| ) (without the square root, because we work with the squared 
equation), while H^ c ^( q ,p) on the R.H.S. is the order-reduced Hamiltonian of p6| , obtained from the highe r-order 
3PN Hamiltonian derived in 0]. Both sides of Eq. J4.25D are written in terms of g and p' by using Eq. ( |4.24| ) . This 
procedure has been used at the 2PN level in 1 13 1, so that we know the values of a\ = \v, a-i = 0, the 2PN values of the 
metric functions A( q ') and D(q'), as well as the 2PN-accurate generating func tion G(q, p' ) (se e e.g. Eqs. (6.24)- (6.27 ) 
in IP^I ). The new unknowns entering at the 3PN level are: «3, a^, Eq. ( 3.3a ), d^, Eq. ( 3.31 ), Z\, z 2 , z 3 , Eq. ( [3.12 ), 
and the 7 arbitrary coefficients ci, . . . , c-j entering the generic form of G3PN: 



G 3 PN( q ,p) = (q- p') 



ci 



p' 6 + -(c 2 p' 4 + c 3 P ' 2 (n • p') 2 + c 4 (n • p') 4 ) + \(c 5 p' 2 + c 6 (n • p') 2 ) + ^ 



g q* q' 

The basic input for writing these equations is the explicit value of the 3PN-accurate Hamiltonian []l4l , ^6|Jl5| 

#rcai(q, p) = #N(q, p) + ffiPN(q, p) + #2PN(q, p) + #3PN(q, p) , 

where 



(4.26) 



(4.27) 



#n (q, p) 



q 



i?iPN (q, P ) = g(3i/ - l)(p 2 ) 2 - \ [(3 + ^)p 2 + Kn ■ P) 2 ] £ + ^5, 



(4.28a) 
(4.28b) 



i?2PN (q, p) = ^ (1 - 51/ + 5^ 2 ) (p 2 ) 3 + J [(5 - 2(V - 3^ 2 ) (p 2 ) 2 - 2^ 2 (n • p) 2 p 2 - 3^ 2 (n • p) 4 ] l - 
16 8 g 



~ [(5 + 8^)p 2 + 3i/(n • p) 2 ] -i - i(l + 3i/)i, 
2 g^ 4 g J 



(4.28c) 



# 3 pn (q, P) = ^ (-5 + 35^ - 70^ 2 + 35^ 3 ) (p 2 ) 4 

+ — \(-7 + 42v - 53i/ 2 - 5i/ 3 ) (p 2 ) 3 + (2 - 3t/W 2 (n ■ p) 2 (p 2 ) 2 + 3(1 - z/V 2 (n ■ p) 4 p 2 - 5!/ 3 (n ■ p) 6 l - 
16 g 



— (-27 + 136^ + 109^ 2 ) (p 2 ) 2 + -^(17 + 30z/Wn ■ p) 2 p 2 + i(5 + 43i/)i/(n • p) 4 
16 lb 12 



25 
'~8~ 



64 



335 
~48~ 



23 



85 



7 



7T z/(n • p) > -r 

16 64 4 I ^ ( „ a 



1 



109 21 2 

-jj _ 32^ + Static j '' 



(4.28d) 



As explained in Refs. {20 16 Tf| and at the beginning of the present section the 3PN Hamiltonian contains one 
dimensionless ambiguity parameter w s tatic- 

When written explicitly, the constraint equation ( 4.25| ) (truncated at 3PN accuracy) yields a system of 11 equations 
for the 10 + 3 unknowns (0:3, 04, (I3, c\, . . . , C7; Z\, z%, Z3) (astatic being assumed to be known). This system can be 
decomposed into three subsystems. The first subsystem consists of 5 equations: 



a 3 + 16ci 
a 3 + 2ci - 2c 2 



-v - iv 2 - 5i/ 3 , 



17 19 2 27 , 
6ci + c 2 - 3c 3 = — v + —v + —v , 
16 4 16 



(4.29a) 
(4.29b) 

(4.29c) 
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o r 3 27 2 81 3 
3C 3 — OCA = v v V , 

2 2 16 ' 



3 9 To 

c 4 = -v* H v 6 . 

2 16 



The second subsystem contains 3 equations: 



3 9 19 

3a 3 + 2c 2 - 2c 5 + zi = -v - -v + — v , 



61 



11 



11 



8c 2 + 6c 3 + 4c 5 — 6c 6 + z 2 = —v + —v — v , 



79 

4c 3 + 10c 4 + 8c 6 + z 3 = ——v 

6 



55 



39 



and the third one consists also of 3 equations: 

2«3 + c 5 - C7 = 

-d.3 + 4c 5 + 6c 6 + 6c 7 = 

— 2a 3 + a 4 + 2c 7 = 



/271 



V 48 + 64* * ]V+ ^ 



5 o 5 o 



35 _3_ 

~8 ~ 32 



7T 2 11/ 



57 



V + 2u A , 



/221 21 2 \ 3 2 1 

-T7T - 777^ + 2 ^static }V+-V + - 



V 12 16 



(4.29d) 
(4.29e) 

(4.30a) 
(4.30b) 
(4.30c) 

(4.31a) 
(4.31b) 
(4.31c) 



The first subsystem, Eqs. (4.29), yields 5 linear equations for the 5 unknowns c\, C2, c 3 , c 4 , and a 3 . It is easily 
found to have a unique solution, namely 



a 3 = 0, 



ci = - — (1 + 3i/ + 51/ 2 ) 1/, 



C2 



C3 



1 

16 



(l + 2z/- lliP)v, 



- — (12 + 48^ + 23^ 2 ) 1/, 
24 v y 



c 4 = — (24 + 7^) v 2 . 
16 



(4.32a) 
(4.32b) 

(4.32c) 

(4.32d) 

(4.32e) 



As already me ntion ed above, note the remarkably simple result a 3 = (which confirms that the energy map takes 
the nice form ( 3.14 )). It is also remarkable that the result a 3 = holds independently of any assumption about the 
"quartic" parameters z\, Z i, and z 3 . 

The second subsystem (4.30) can be viewed (after inserting the unique solution of the first subsystem) as an 
overdetermined system for the two unknowns C5, eg. It is then easily seen that it will admit a solution if and only if 
the parameters z 1; z 2j and £ 3 satisfy the following linear constraint: 



8zi +4z 2 + 3z 3 = 6(4-3^)2^. 



(4.33) 



This linear constraint forbids us to consider the simplest "geodesic" case where z\ — z^ = Z3 = 0. We can, however, 
continue to impose the natural conditions z\ = — zo which simplify very much the 3PN effective dynamics of circular 
orbits. With this choice, the general constraint ( 4.33 ) yields 



z 3 = 2(4 - Zv)v. 

Having fixed the values of z±, z 2 , and z 3 , the system ( 4.30| ) uniquely determines C5 and cq to be 



(4.34) 



16 



c 5 = - — (13-16i/ + 6v 2 )v, 



c 6 = (115 + 116f - 26i/ 2 ) i/. 



(4.35a) 
(4.35b) 



Finally, the subsystem ( |4.3l| ) gives 3 equations for the remaining 3 unknowns cj, g?3, and a 4 . The unique solution 
of this system reads: 



C7 



1 , 155\ 3 2 1 o 

64^ + *rr + r + r < 



d 3 = 2(3f - 26) ^, 



«4 



94 41^ 
3~ ~ 32' 



— — 7T 2 + 2(jJ sta tic I ■"■ 



(4.36a) 



(4.36b) 



(4.36c) 



No te how si mple the structure of the coefficient a 4 , Eq. (4.36c), is. Indeed, the right-hand-sides of all the equations 
( 4.29 ), ( 4.30 ), ( 4.31 ), were polynomials of the third degree in v. Therefore, one would have a priori expected the 3PN 
coefficient a 4 to have the same structure: 04 = an v + 042 v 2 -\-a/^v i . It is remarkable that the coefficients of v 2 and v 3 
happen to vanish in 04 (such a simplification does not occur in the 3PN-level coefficients appearing in the functions 
E(x), e(x) and j 2 (x) considered above). This simple structure of 04 can be brought out by defining the following 
quantity, 



3 ' 64 

in terms of which the value of 0,4 can be written as 



-? + ^ 2 = -9.3439. 



a 4 = 2(u> s 



(4.37) 



(4.38) 



The presence of already two cancellations in 04 

Eq. d4.36c[) 



so that 



042 = = 043) suggests that the yet undetermined value of w s tatic 
simply vanishes. We shall see below that this conjecture is 



might be precisely ui* t 

indirectly supported by the fact that a numerical value w s tatic — — 9 is selected by the requirement that the various 
methods discussed in this paper agree in their predictions for LSO quantities. Note also the remarkable fact that if 
^static = ^static au tne 71-2 t erms cancell in all the 3PN-level coefficients, so that they all become rational (as were the 
2PN ones). [Stated in reverse, if one could a priori prove that all the 3PN coefficients are rational, this would support 
the conjecture that w s t a tic = Static > though it would be also compatible with having w s t a tic = w static+ a rational 
number.] 

The coefficient 04 enters the PN expansion of A(u) = —.900(9') (with u = 1/V): 



A(u) = 1 - 2u + 2vv? + a 4 (v)u 4 + 0(u 5 ). 



(4.39) 



As mentioned above, we improve the behaviour of the PN expansion of A(u) by Padeeing it: A Pn {u) = Pf[T n+l [A{u)]] 
with k + £ = n + 1. We impose the constraint k > to inject the information that A(u) should qualitatively look like 
j4s c hw( w ) — 1 — 2m, i.e. have a zero at some u = | + 0{v). As said above, the most robust (for our purpose) Pades of 
A(u) are the ones with k = 1 and £ = n. Finally, we get the following sequence of Pade-improved A: 



A Pl (u) = Pl[T 2 [A{u)]] = l-2u, 



A P2 (u) = Pi [T 3 [A(u)}} = 



l-(2 



1 + -zVU-\- vu 2 



A P3 (u) ee Pi [T 4 [A(u)}} = 



2(4- v) + (04(f) - 16 + 8i/) i 



2(4 - v) + (a 4 (u) + Av)u + 2(a 4 (v) + Av)u 2 + 4(a 4 (i/) + v 2 )u 3 ' 
To extract the LSO quantities from these Padeed A's, we must consider the effective radial potential 



(4.40a) 
(4.40b) 

(4.40c) 



W^(u)=A P Ju)(l+j 2 u 2 ) 



(4.41) 
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The value of j for which this radial pote ntial h as an infle ction point defines j'lso { v ) \ the corresponding value of u being 
wlso = u *0lso)- As explained in Eqs. ( 3.21 ) and ( 3.22 ) above one then deduces the energy and the orbital frequency 
of the LSO. Note that the 2PN Padeed Ap 2 that we use here differs from the straightforward Taylor approximant At 2 
used in Ref. p3| . However, this difference is essentially negligible (as shown by comparing the lines "eff. method" 
and "BD" in Table I). The Pade improvement is, however, rather important at 3PN in the case where w s t a tic is 
significantly larger than w* tatic . Indeed, in this case, Padeeing allows one to tame the effect of a largish (positive) 
0,4(1;): 0,4 (4j ~ 4.67+ ^static- For instance, when v = j and w s t a tic ^ —1.2 the radial potential built from a straight 
Taylor-approximated (u) would not give rise to an inflection point continuously connected to the test-mass limit. 
We attribute this lack of structural stability to the known bad properties of high-order PN expansions, and not to 
the effective-one-body approach. The (numerical) results obtained by the effective one-body approach are exhibited 
in Table § 



V. DISCUSSION 



Before discussing the meaning of the results obtained above, let us state what we would a priori expect. First, we 
recall that the study in Ref. JlJ has shown that the sequence of Pade approximants of the invariant function F(v), giving 
the gravitational wave flux in terms of v = (GMlu/c 3 ) 1 ^ 3 = x 1 / 2 , had very good (and very monotonic) convergence 
properties toward the exact result. [By contrast, the sequence of Taylor approximants was badly convergent, and 
unstable when v < VlscT = 0.40825; see Figs. 3a and 3b of jjj.] In our case, as one can meaningfully (at least for 
the j- and effective-one-body methods) consider the 1PN, 2PN, and 3PN approximations, we would expect that a 
good resummation technique would ensure that any LSO quantity, say QlsOj be determined with increasing accuracy, 
when using higher PN information, and, more precisely, that 

QfSo^LSO + a(^L S or + \ (5.1) 

with (hopefully) coefficients a and b small enough to ensure a visible convergence (when xlso — ^lscT — 1/6)- As a 
minimum test of improved convergence we hope that IQ3PN — <32Pn| would be significantly smaller than IQ2PN — Qipn|, 
i.e. that the addition of the 3PN information would have only slightly refined the previous 2PN estimates of LSO 
quantities [p],^3|. 

Independently of this expectation, we had also hoped, when starting this investigation, that the LSO quantities 
might be "robust" under the lack of precise knowledge of a sole ambiguous coefficient (astatic) among many others 
in -/?3pn- [Given that the amplitude of this coefficient would have some plausible upper bound; as discussed in the 
Appendix A of [jig] .] However, the results exhibited in Table | and Fig. [jjshow that, in spite of our use of resummation 
techniques, the LSO quantities appear to be quite sensitive to the exact value of w s tatic- A first conclusion of our 
work is therefore that it is quite important to resolve the problem of static ambiguity, arising at 3PN when using 
delta functions to represent compact (but extended) objects. Until this problem is unambiguously solved, it will not 
be clear whether (as proposed in [[[]|l|) it is possible to trust suitably resummed versions of PN-expanded results. 

In the meantime, however, we wish to point out several remarkable features of the dependence of our various results 
on ^static- In Fig. |2] we plot (for the equal-mass case, v = 1/4) our various predictions, at the 3PN level, and using 
various methods, as a function of the 3PN ambiguity parameter ^static- It is quite interesting to note that two a priori 
independent things happen: 

(i) there is a value of cj sta tic , namely 

Otitic * "9, (5.2) 

for which the three different methods give, at 3PN, nearly coincident LSO predictions. 

(ii) For this "best fit" value tj^^ ic the 3PN LSO predictions exhibit the expected convergence property that 
IQ3PN — Q2Pn| is significantly smaller than IQ2PN — Qipn| (see Table H below). 

We have checked that these two remarkable properties hol d for all values of the parameter v < j. Actually there 



are several other ways of selecting the appro xima te value (5.2), i.e. of understanding why it plays a special role. First, 



we have seen above that the precise value ( 4.37 ) (which is near the "best fit" value (5.2)), played a special role in 
simplifying not only 04 but also all the other 3PN coefficients. Second, it seems natural to expect that the true value 
of (^static will be such that the full Taylor expansions of most of the invariant functions will be smooth deformations 
of their test-mass limits. A minimum requirement for this property of "structural stability" under the turning-on 
of the parameter v seems to be that the Taylor coefficients of the functions e(x), j 2 (x), K(y), and 1/A(u) do not 
change sign as v varies from to 1/4 (we only consider functions with infinitely many non-zero Taylor coefficients in 
the test-mass limit). [One could actually impose a more restricted bound on the (/—variation of the 3PN coefficients, 
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FIG. 2. Equal-mass (y — 4) binary systems LSO parameters obtai ned by means of different methods as functions of the 
ambig uity parameter aj sta ti c - The reduced binding energy -Elso, Eq. (2.2), and the reduced angular momentum j'lso, Eq. 
( 2.18 ), are divided by their "Schwarzschild values": |i?LSo'| = ^ ~ ant ^ jtso — V^12; the dimensionless orbital frequency 
lolso is defined in Eq. (|f.9|). The lines shown in the plots correspond to different methods: j'-method ( solid) , effective one-body 
method (dotted), and e- method (dashed). The solid vertical lines correspond to w s t a tlc = ^static i Eq. (4.37). 



especially given the information that the 2PN coefficients are found to vary by a smallish fractional amount.] We find 
that this minimum requirement is satisfied only if w s t a tic < —0.62 (the consideration of the expansion of e(x) gives the 
most stringent bound). Another natural requirement for "structural stability" under v ^ would be to impose that 
all the near-diagonal Pades (without restricting oneself, as above, to the most "robust" ones) of the functions above 
exhibit real poles that are smoothly connected to their test-mass counterparts. This requirement is most strin gen t 



when considering P%(\/ A(u)), and yields the limit astatic < —8.35. Combining this with the general limits Eq. (4.1) 
suggests that w s t a tic lies within the small range —10 < ^static < —8.35 



Let us quote a last way of selecting the value (5.2). It consists in comparing the 3PN Taylor coefficients of the 
invariant functions that contain only a finite number of terms in the test-mass limit. For instance, consider A(u) and 
l/j 2 (x). In the test-mass limit A(u) = 1 — 2u and l/j 2 (x) — x — 3x 2 . When v ^ there will appear further powers of 
u or x with coefficients vanishing with v. At the 3PN level there is a term a^u A in A(u), and a term |&4X 4 in l/j 2 (x). 
[The factor 8/3 is introduced to have the same (linear) dependence on w s t a tic in «4 and 64.] These two coefficients 
read 

/94 41 \ 

04 (f) = ( -3- _ l^ 11 J v + 2 ^static V 

w 18.6879 v + 2 u; sta tic v (5.3) 

and 



, . , /5269 41 2 \ 61 2 1 , „ 
y ' \ 192 32 J 32 216 



14.7973 v - 1.9063 v z + 0.00463 v A + 2 w sta tic v. (5.4) 



Let us first note that the terms oc v 2 and v 3 in Eq. (5.4) are numerically nearly negligible. Forgetting about them 



(i.e. working with b' 4 = (^Pr _ % 1T ' 2 ) v )i we t nen see, by comparing Eqs. (|5.3[) and (^3.4[), that the coefficients an{v) 



and b'^v) are approximately identical. In particular, this means that there will be a small range of values of w stat j c for 
which ci4 and b\ will be simultaneously small. The existence of this range explains why the j- and effective one-body 
methods can give numerically similar results. We can then make an analytical estimate of the 'best' value of the 
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TABLE II. Equal-mass (y = 1/4) binary systems LSO parameters obtained by means of t he j- and effective-one-body 
methods. The 3PN values of the LSO parameters were calculated for astatic = w* tatic , Eq. (4.37); Eq. (5.5) defines the 
S-estimates. 



method 


IPN 


2PN 


/I ijiSchw 1 

/l-^LSO 1 

3PN 


S-estimate 


IPN 


/,-Schw 

Jlso/Jlso 
2PN 3PN 


S-estimate 


IPN 


2PN 


WLSO 

3PN 


S-estimate 


j-method 
eff. method 


-0.973 
-1.007 


-1.091 
-1.048 


-1.054 
-1.049 


-1.063 
-1.049 


1.014 
1.000 


0.970 0.982 
0.983 0.983 


0.979 
0.983 


0.960 
1.015 


1.173 
1.075 


1.091 
1.077 


1.114 
1.077 



ambiguity parameter w s tatic by looking for the value of w s tatic which simultaneously minimizes (in a least square sense) 
ct4 and 64. It is easily seen that the expression [ci4(^)] 2 + [fr 4 (z^)] 2 attains its minimal value (as function of astatic) for 



41 2 11285 
64 71 " ~ 768 



- -< - ' - ~ ^rF ~ " 8 - 37 • ( 5 - 5 ) 



This numerical value is not very far the special values selected by the other arguments discussed above. 

Summarizing: several (partially) independent arguments suggest that the true value of astatic lies in the range 
— 10 < (^static ^ — 8. For definiteness, and for the purpose of the following discussion, we shall henceforth assume that 
the "correct" value of w static is 

^static = Static- (5-6) 

In Fig. ^| we have included vertical lines corresponding to w s tatic = w* tatic , to show visually that Eq. ( |5.6| ) is well 
compatible with our argument based on the convergence of the various methods. 

One can also see in Fig. || that the curves related to the e- and j-methods have a second intersection point, besides 
the one around w^atfc — — 9- However, we have checked that the value of astatic at this point strongly depends on the 
value of the parameter v. For this reason, and also for the fact that this point does not give an agreement with the 
"effective" method, we do not take this second intersection point as evidence for a different value of w s tatic 

Admitting (for the sake of the following argument) Eq. ( |5.6| ) we wish to propose a further way of improving the 
accuracy of the predictions of LSO observables. Indeed, if one has at one's disposal three successive approximations, 
namely the IPN, 2PN, and 3PN estimates of some quantity Qlso, one can combine this information to refine the 
estimate of the (unknown) exact value Qlso- The rationale for this is to assume that the approach to the limit, when 



the order n of the approximant increases, is approximately described by Eq. (5.1) (i.e. essentially that the inaccuracy 
of the nth estimate decreases proportionally to the (n + l)th power of some constant c = 6xlso < !)• Then, under 
this assumption the knowledge of three (successive) approximants, say Q n -i, Q n , and Q n +i, gives three equations 
(Qm = Qx + ac m+1 ) for the three unknowns (Qx,a,c). One can solve this system of equations and deduce, in 
particular, the value of the looked for n — > oo limit Qx in terms of Q n -\, Q n , and Q n +i- The result defines the 
so-called "Shanks transformation" |2^], namely 

Qx ~ S [Q] = On— l ~ Qn j-g 

Qn+l + Qn-1 ~ 2 Q n 

When one disposes of more than three Q m 's, the Shanks transformation associates to the original (truncated) sequence 
{Qii Q2, ■ ■ ■ , Qn) a shorter, but hopefully faster converging sequence (S2 [Q], ■ ■ ■ , Sn-i [Q])- In our case, the Shanks 
procedure associates to any triplet of LSO quantities (Qi, Q2, Q3) = (Qipn; Q2PN' Q3PN) a single number, 

r>Lso r>Lso (n^so\2 
— n LSo 1 rjLSo _ 9 (qlso ' y°-°> 

^3PN "T" ^1PN z ^2PN 



which is a (hopefully better) estimate of the (unknown) exact value Qlso- We shall refer to (5.8) as the S'-estimate 
ofQ LS °. 

In Table [jl] ( see also Fig. |) we ap ply this procedure to our two best methods: the j-method and the effective- 
one-body one, under the assumption ( |5.6| ) (which is needed to exhibit a visible convergence among the first three PN 
approximations) . 

Given our present (incomplete) knowledge we consider that the S'-estimates exhibited in Table || represent our 
best estimates of LSO observables. To verify the plausiblity of these estimates one should resolve the issue of the 
ambiguous coefficient w s t a tic in the 3PN dynamics. [In principle, this can be done by implementing the matching 
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FIG. 3. Equal-mass (y — 1/4) binary systems reduced binding energy i^Lse/I^Eso"! versus the dimensionless orbital 
frequency o)lso, f° r different methods discussed in our paper. For j- and effective one-body methods we have plotted the 
results at the 1PN level, 2PN level, and S-approximants; they are all exhibited in Table [n|. We have also shown the results 
obtained in Refs. [8] (labelled by KWW), and by applying the 2PN effective-one-body method to the "Wilson-Mathews" 
truncation of general relativity (labelled by WM, see Appendix B). 



meth od described in [£3| and used there at the 2PN level.] If this resolution approximately confirms the estimate 
( |5.6|) the 5-estimates will be confirmed. If a very different value of ast atic i s found, it might still be compatible with a 
less evidently convergent PN sequence. And hopefully, the S'-estimate (5.8) of this new sequence will give an improved 
3PN-accurate estimate of LSO observables. 

Under the assumption that the iS-estimates are accurate, there are several interesting conclusions that we can 
draw. First, we remark that the final estimates are quite near the 2PN- level predictions of the effective one-body 
approach, see Fig. |^. Although this may seem disappointing (an enormous, not yet completed, 3PN work leading to 
a confirmation of 2PN estimates), this would be a scientifically very useful conclusion. Indeed, this would (in our 
minds at least) establish the soundness of the philosophy advocated in Refs. |l],[l3| and here, namely that resummation 
methods can be meaningfully employed to make analytical predictions concerning physics near the Last Stable Orbit. 
This would then also give support to the recent work of Buonanno and Damour |24|] in which the 2PN effective 
one-body Hamiltonian has been used, together with Pade-resummed estimates of gravitational-radiation damping, to 
study the transition between the inspiral motion and the final plunge of a binary system. [Let us note, in passing, 
that this work shows that, though it is crucial to have good initial estimates of the LSO quantities defined by the 
Hamiltonian, the final observable effects linked to the presence of an LSO are blurred by radiation-reaction effects.] 
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APPENDIX A: 3PN EFFECTIVE "GEODESIC" ONE-BODY DYNAMICS 

We consider here an effective 'relativistic' one-body Hamiltonian of the simple "geodesic" form 
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H*ti,p') = JA( q >) 



i+p,2+ (^H (n '- p,)2 



(Al) 



The Hamiltonian is related to the real 'non-relativistic' Hamiltonian H™£, Eq. (4.27), through the constraint 
equation 



H*( q >( q ,p'),p')] 2 = [l + H™(q,p(q,p>)) [l + ai + a 2 + a 3 (jf£2) 



(A2) 



As in the text both sides of Eq. (A2) are written in terms of the variables q and p' by means of Eqs. (4.24) with a 
generating function G of the form 



where 



G(q,p') = Gip N (<7,p') + G 2 pn(9,p') + G 3PN (q,p'), 



Gn-Nl/yV) M l' (,9iP' 2 + | 



G 2 pn(<?,p') = (q • p') 



G 3 PN(q,p') = (q • p') 



b lP '* + -(b 2 p ,2 + 6 3 (n •p') 2 ) + ^ 



ci p' 6 + - (c 2 p' 4 + c 3 p' 2 (n • p') 2 + c 4 (n • p') 4 ) + 4(c 5 p' 2 + c 6 (n • p') 2 ) + % 



(A3) 

(A4a) 
(A4b) 
(A4c) 



Written explicitly, the constraint equation (A2) is equivalent to a sytem of 23 equations for the 23 unknowns: oi, 
. . , a^] di, d2, d 3 \ ai, a 2l a 3 ] gi, g 2 ] b\, . . . , 64; C\, . . . , C7. The unique solution to these equations reads 



01 = -2, 



«2 



a 3 



04 



(3i/ -4)i/ 
2(5 -2i/)' 

(1600 - 1576i/ + 392j/ 2 - 9i/ 3 )i/ 
16(5- 2i/) 2 ' 

(4280 - 3349^ + 692j/ 2 - 9^ 3 )zy 41tt 2 



6(5-2^) 2 
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V + 2w sta tic V, 



ai 



a 2 



a 3 



(3f - 4)v 
2(5 - 2u) ' 

(-2400 + 1936t/ - 408t/ 2 + 9^ 3 )^ 
16(5- 2v) 2 ' 

(-486400 + 703680j/ - 383904^ + 93704i/ 3 - 8580i/ 4 - 27i^ 5 )^ 

64(5- 2vf ' 

(Av - 2>)v 
2(5 -2vY 

(80- 32?/ + 7^ 2 )(3t/ - A)v 
4(5 - 2i/) 2 ' 

(3650 - 3660i/ + 1829j/ 2 - A21v 3 + 44z/)(3z/ - A)v 
8(5 - 2^) 3 ' 



(A5a) 
(A5b) 

(A5c) 

(A5d) 

(A5e) 

(A5f) 

(A5g) 

(A5h) 

(A5i) 

(A5j) 
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(. - 6). 

51 = W^ty (A5k) 

5(4 - 2. + . 2 ) 

92 = 4(5-2.) ' (A51) 



(210 - 178. + 50. 2 -3. 3 ). 
16(5- 2.) 2 



5 i = — — — > ( A5m ) 



(-350 + 419. - 149. 2 + 13. 3 ). /4 . 

62 = " 8(5^2^ (A5n) 



(120 -81. + 38. 2 -6. 3 ). 



8(5-2.) 



2 



200 + 200. - 816. 2 + 360. 3 - 49. 4 
32(5 - 2.) 



(A5o) 



°4 = 5^ ?TA2 ' ( A5 P) 



(5900 - 12700. + 8270. 2 - 2082. 3 + 192. 4 - ll. 5 ). . . „ , 

Cl = 64(5^ ' (A5q) 

(-25850 + 44320. - 24876. 2 + 5157. 3 - 283. 4 + 13. 5 ). . . . 

C * = - 32(5^ -' (A5r) 

(-25200 + 21360.- 4364. 2 - 1077. 3 + 408. 4 - 47. 5 ). . 

C3 = 96(5^ ' (A5S) 

(2160 - 1834. + 900. 2 - 228. 3 + 23. 4 ). 2 x 

C4 = ' 32(5^^ (A5t) 

(247000 - 407080. + 225416. 2 - 49240. 3 + 3730. 4 - 151. 5 ). , , . 

05 = 128(5 -2.) 3 ' (A5U) 

7(10300 - 16360. + 7136. 2 - 1372. 3 + 154. 4 - 7. 5 ). 

C6 = 192(5 -2.) 3 < (A5V) 

(-962800+ 1472880. - 813024. 2 + 206500. 3 - 25110. 4 + 1479. 5 ). it 2 N 

C7 = 384(5 - 2.) 3 64" (A5w) 

As said in the text, in view of the complexity of these results, we do not take this possibility seriously. We prefer to 
it the non-minimal ("non-geodesic") Hamiltonian given in Sec. IV. 



APPENDIX B: 2PN RESULTS FOR THE CONFORM ALLY-FLAT TRUNCATION OF GENERAL 

RELATIVITY 



By contrast, let us note that other approximation philosophies are, in our opinion, less reliable to make predictions 
concerning the LSO. We have in mind here: (i) the use of non-resummed (or only partially resummed) PN expansions, 
and (ii) the " Wilson-Mathews" -type |2f| truncation of Einstein's theory, in which the spatial metric is taken to be 
conformally flat. As an example of the first philosophy, let us consider the proposal of Kidder, Will, and Wiseman |ic| l 
to partially resum the Damour-Deruelle equations of motion by separating out (and resumming) the "Schwarzschild" 
(. = 0) terms. This approach led to the prediction (at 2PN) that the LSO is significantly less bound (when . = 1/4) 
than the "Schwarzschild" limit. In terms of orbital frequency at the LSO, Ref. |l(J predicts u3lso(1/4) — 0.891 < 1. 
This contrasts very much with our 2PN and 3PN estimates above which consistently indicate that wlso(^) is larger 
than one (and that the LSO is more bound than its Schwarzschild limit: ^lso/I^lscTI < — -0- Independently of 
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this (biassed) argument, we think that both the robustness and the consistency of the "hybrid" approach of [[l0| are 
seriously in doubt. Indeed, Refs. 11 and fl^] have shown that the hybrid approach was robust neither under a change 
of formulation (Hamiltonian versus equations-of- motion), nor under a change of coordinate system. Moreover, Ref. 
jjj has questioned the consistency of this approach by pointing out that the non-resummed V-corrections" represent, 
in several cases, a very large (larger than 100%) modification of the corresponding ^-independent terms. 

Regarding the conformally-flat truncation it was noted by Rieth |^6| that this implies significant deviations from 
the Einstein dynamics already at the 2PN level. We have investigated this question further. In Ref. we gave the 
invariant functions defined (at 2PN accuracy) by the Wilson-Mathews-type truncation. Applying now our j-method, 
we find (at 2PN) 



1 



1 + I (9 + v) x + -!- (324 - 333i/ - 50i/ 2 ) x 1 
3 So 



(Bl) 



whose Padccd form is 



JWMP 2 ( X ) = P l [^[j'wM^)]] = 



1 + 5"+(t* 



'(1 



37, 



(B2) 



This leads to a prediction for the 2PN £lso which can be written down analytically 



ox LSO 



8(9- 



3(15 + 2i/)i/ 



2(9 



V324 - 333z/ - 50i/ 2 



(B3) 



The corresponding dimensionless orbital frequency, for v = 1/4, equals ljlso = 1-4378, the reduced binding energy 
£lso/|-Elso w I = -1-2253, and the reduced angular momentum jLso/jEscT = 0.9293. 

We have also studied, at the 2PN level, the effective one-body method for the Wilson-Mathews dynamics. Using 
the procedure decribed in Sec. IV D above, imposing at the 1PN level the condition di = 0, we have found that the 
effective-metric function Awm(u) at the 2PN accuracy reads 



A wu (u) = l-2u+ a 3 {is)u 3 + C(m 4 ), 



(B4) 



where 



a 3 (v) = -(18-5i/)i/. 
We have improved the behaviour of the 2PN expansion of Avvm(u) by Padeeing it: 



(B5) 



lWMP 2 



(u) = Pl [T 3 [A WM (u)}} = 



16 



1 



16' 



5„2\ 



(B6) 



To extract the LSO quantities from this Padeed A, we have considered the inflection point of the effective radial 
potential Awmp, + 3 u \ w hich defines the angular momentum j'lsoC^) and the location wlso(i/) of the LSO; 
then using Eqs. ( 3.21 ) and ( 3.22 ) one calculates the energy and the orbital frequency of the LSO. The results, for 
v = 1/4, are: dimensionless orbital frequency Qlso — 1-1482, reduced binding energy ^lso/I^lscTI = — 1-0972, and 
reduced angular momentum Jlso/JlscT = 0.9647. Let us also mention that Cook, using another conformally flat 
approximation 0, obtained, for v = 1/4, the following LSO parameters: the dimensionless orbital frequency wlso = 

1.579, and the reduced angular momentum j'lso/j'lscT ~ 0.8591. 
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